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A Land Climate Data Record 

Multi instrument/Multi sensor Science Quality Data Records used to 
quantify trends and changes 
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Emphasis on data consistency - characterization 
rather than degrading/smoothing the data 
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Land Climate Data Record (Approach) 


Needs to address geolocation, calibration, atmospheric/ BRDF correction issues 
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NDVI at California Redwood Site. 1951-1999 
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MODIS used as a reference for past and future 
land data record (example NDVI) 


Evaluation overAERONET (2003) 
0.007 <Precision < 0.017 
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Independent evaluation of the precision 
Over 2000-2004 CMG daily time series 
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MODIS product and validation methodology 
used to evaluate other surface reflectance 
product: example LANDSAT TM/ETM+ 


• WELD (D. Roy) 120 acquisitions over 23 AERONET 
sites (CONUS) 

Junchang Ju, David P. Roy, Eric Vermote, Jeffrey Masek, Valeriy Kovalskyy, Continental-scale validation of 
MODIS-based and LEDAPS Landsat ETM+ atmospheric correction methods, Remote Sensing of Environment 
(2012), Available online 10 February 2012, ISSN 0034-4257, 10.1016/j.rse.201 1.12.025. 


• GFCC: Comparison with MODIS SR products 

- GLS 2000 demonstration 

Min Feng, Chengquan Huang, Saurabh Channan, Eric F. Vermote, Jeffrey G. Masek, John R. Townshend, 
Quality assessment of Landsat surface reflectance products using MODIS data, Computers & Geosciences, 
Volume 38, Issue 1, January 2012, Pages 9-22, ISSN 0098-3004, 10.1016. 

- GLS 2005 (TM and ETM+) 

Min Feng Joseph O. Sexton, Chengquan Huang, Jeffrey G. Masek, Eric F. Vermote, Feng Gao, Raghuram 
Narasimhan, Saurabh Channan, Robert E. Wolfe, John R. Townshend .Global, long-term surface reflectance 
records from Landsat: a comparison of the Global Land Survey and MODIS surface reflectance datasets. 

Remote Sensing of the Environment (in review) 


Land Product Validation and Evolution Workshop, Frascati, Italy, January 28-30, 2014 


WELD/LEDAPS results (Red-band3) 

LEDAPS WELD uses MODIS aerosol 




Top of the atmosphere 



Surface Reflectance Truth 
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Visible Near IR 


MODIS Reflectance time series 


• Reflectance time series show high-frequency variability 

• The “noise” is partly due to directional effects. 

• Selection of specific geometries decreases temporal coverage 

• Can we correct for the directional effect and retain the original 
temporal resolution ? 
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Analytical model and correction 

Linear models : 

P (0 S > 0 V , ♦) = K + k i F i(0 s , B v , <|>) + k 2 F 2 (0 s , 0 v , <|>) 

Several choices for Fi and F 2 

Fj: Model surface effects (soil roughness) 

F 2 : Model volume effects (R.T. within canopy) 

p (0 S , 0 V , ® = k 0 (t) [ 1+ kj/ko Fj(0 s , e v , ® + k 2 /k 0 F 2 (0 s , 0 v , <|>) ] 

= k 0 (t) [ 1+ R F,(0 S , 0 V , <►) + V F 2 (0 s , 0 v , <|>) ] 


Correction: 

P cor = p(0 8 , 0 V , 4>) [ 1+ R F 1 (45,0,0) + V F 2 (45,0,0) ] / 

[ 1+ R Fi(6 s , 0 V , 4>) + V F 2 (0 s , 0„ <(>) ] 

Maignan et al., Rem. Sens. Env., 2004 



Which is the “best” model ? 




Cumulative histogram of reflectance [%] error of fit. 

• Look for the parameters for a best fit (invert the model) 

• Compute error of fit 

• Among the 6 tested models, RossLiHS allows the best fit 

• Clear improvement when using Hot Spot correction 


Maignan et al., Rem. Sens. Env., 2004 


BRDF model inversion 
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R and V are linear function of the NDVI 

We look at the difference between 
successive measurements 

Notations used here for an easy 
inversion of the model parameters 

Matrix writing: 
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Data location 
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• MODIS data are distributed as “tiles” (10° of lat.) 

• To limit data volume, we focus on a single tile 

• Select a tile over Eastern Australia for (i) variety of surface 
cover, (ii) number of clear observations, (iii) low aerosol load 



R parameter Channel 2 


BRDF parameters: R and V 
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Visible band 
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Near IR band 


Analyzed the BRDF parameters distributed in the official 
MODIS products 

Parameters show very unrealistic temporal variations. 
Our method shows more realistic results 


Quantification of time series noise 

• For each triplet of observations, one can estimate middle one 
from the earlier and later: 

* (ti ti-i)Pi+i "I" (tj+i ti) Pi— l 

Pi — 7 _ 

L l+1 L l-1 


One can then compute a “noise” from the quadratic sum of the 
difference between the measurement and their interpolated 
counterpart: 


ixy t (Pt- Pi y 

<j 2 (p) = — i+i i ~\ 

yN-l ± 

Ll = 2 t i + 1 - t(_! 

We use this definition in the following to quantify the time series 
quality 


5x5 resolution Noise 5x5 resolution Noise 



Impact of spatial scale 




1x1 resolution Noise 1x1 resolution Noise 


•The noise of the corrected 
time series is much larger than 
that we obtained earlier using 
CMG (Climate Modeling Grid : 
5 km) lower resolution data. 

•We show here a comparison 
of the noise obtained at the full 
resolution against that 
obtained when aggregating 
5x5 pixels. 
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Noise vs Spatial heterogeneity 
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•There is a very strong 
correlation between the 
spatial heterogeneity 
(quantified here as the 3x3 
standard deviation) and the 
noise on the corrected time 
series. 

•Clearly, the spatial 
heterogeneity affects the 
quality of the time series and 
there is an easy explanation 
for that 
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Impact of spatial scale 



The “noise” of the time 
series decreases when 
the spatial aggregation 
increases. There 
seems to be an optimal 
scale at 2 km (4x4 
pixels) 


Conclusions (1/2) 

• Directional effects on the Earth reflectances are large (factor of 
2 to 4 depending on wavelength) 

• There are simple analytical models (3 parameters) that 
reproduce accurately the observed signatures 

• The reflectance is modelled as the product of a normalized 
reflectance, that may vary rapidly, and a BRDF model (2 
parameters) that varies more slowly. 

• The two model parameters can be parameterized as a linear 
function of the NDVI. 

• We have developed a method to estimate easily these 
parameters from the reflectance time series 

• Corrected time series are much smoother than their original 
counterpart, and can be used to extract fine signal 


Conclusions (2/2) 


• The official MODIS BRDF parameters are unreliable 

• The time series at the full (500 m) resolution appear noisier 
than at lower resolution 

• Spatial heterogeneity of the reflectance is the driving factor for 
this additional noise 

• We suggest an optimal resolution of 2 km for the use of MODIS 
data time series 


Use of BRDF correction for 
product cross-comparison 


Comparison of aggregated FORMOSAT-2 reflectance and MODIS reflectance. No BRDF 
correction. Density function from light grey (minimum) to black (maximum); white = no 
data. 


Comparison of aggregated FORMOSAT-2 reflectance and BRDF corrected MODIS 
reflectance. Corrections were performed with Vermote al. (2009) method using for each day 
of acquisition, the angular configuration of FORMOSAT-2 data. 
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Cross 


Calibration 


of NOAA 16 AVHRR 
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The coefficients 
were consistent within 
less than 1% 
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and band 2 (red line and square). 
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Use of BRDF corrected reflectance for cloud 
mask evaluation of AVHRR Time Series 
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CLAVR (Pathfinder II) cloud mask 
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Use of BRDF corrected reflectance for cloud 
mask evaluation of AVHRR Time Series 


DBQ4 [22.88S, 43.67E] ( 27.35%) 
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Using Direct comparison with MODIS Aqua for validation 


Comparison of MODIS Aqua and 
NOAA16 AVHRR data, A (Red) 
,B (NIR) ,C (NDVI) are observed 
over AERONET sites for 2003- 
2004, D (Red), E(NIR), F(NDVI) 
are simulated using a vegetation 
model that account for spectral 
difference between MODIS and 
AVHRR bands. G shows over 
the AERONET sites MODIS 
NDVI versus corrected AVHRR 
NDVI computed from spectrally 
adjusted AVHRR surface 
reflectance. 
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Simulated data 
y = 0.937 x + 0.006 
RMSD = 0.036 (10%) 
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VI IRS band M7 


Use of BRDF correction (VI IRS) 



VI IRS SR product Aqua SR product 
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VIIRS calibration is being monitored on a 
continuous basis (selected daily obs) 


1*04 
1*02 
1 

0*9B 
0*96 

2012 2012*2 2012*4 2012*6 2012*6 2013 

Land Product Validation and Evolution Workshop, Frascati, Italy, January 28-30, 2014 




VI IRS calibration monitoring — |— 
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Conclusions 


Surface reflectance algorithm is mature and pathway 
toward validation and automated QA is clearly identified. 

Algorithm is generic and tied to documented validated 
radiative transfer code enabling easier inter-comparison 
and fusion of products from different sensors 
(MODIS,VIIRS,AVHRR, LDCM, Landsat, Sentinel 2 ...) 
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